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In this proceedings, I consider two-loop contributions to 
real photon production in thermal QCD. If the photon is 
strictly massless, strong collinear divergences appear in this 
calculation. These singularities are regularized by the quark 
thermal mass of order gT, which generates powers of 1/g, 
so that the corresponding terms are strongly enhanced with 
respect to naive expectations of their order of magnitude. 
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I. INTRODUCTION 

From a phenomenological point of view, photon pro- 
duction is considered as a potential candidate for the de- 
tection and discrimination of a quark gluon plasma. This 
quantity is also interesting from a more formal point of 
view in order to test our control over some aspects of ther- 
mal field theories like infrared and collinear singularities. 
One reason is that this quantity is obviously observable 
and should therefore be infrared and collinear safe in a 
consistent theory. 

At one loop, the thermal production of soft real pho- 
tons |l],|2| suffers from a logarithmic collinear singularity, 
that can be cured by an extension of the hard thermal 
loop resummation technique proposed by Flechsig and 
Rebhan || . It appeared also that common processes like 
bremsstrahlung are not present at one loop, while they 
are assumed to play the dominant role in semi-classical 
approaches 

Therefore, we considered the two-loop diagrams giv- 
ing bremsstrahlung js] |l0| . It turns out that they suffer 
from collinear singularities much stronger than the 1-loop 
ones. After regularization by a thermal mass of order gT, 
these singularities give extra factors 1/g, which spoils the 
natural hierarchy of perturbative expansion. I present 
this calculation in this proceedings. Then, I compare the 
result with that of other approaches, and finally, I con- 
clude with some open questions related to what might 
happen at higher orders. 



*Work done in collaboration with P. Aurenche, R. Kobes, 
E. Petigirard and H. Zaraket. Talk presented at the TFT'98 
conference, Regensburg, Germany, 10-14 august 1998. 



II. TOOLS 

A. Hard thermal loops 

This first section is devoted to a presentation of the 
tools we used in order to perform this calculation. The 
general framework of our study is the effective perturba- 
tive expansion that one obtains after the resummation of 
hard thermal loops Jll],[l^]. To be more definite, we used 
the retarded-advanced version of the real time formalism 
@0. Initially, one of the motivations for hard thermal 
loops was the resolution of the long standing puzzle of the 
gluon damping rate, but they provided also a consistent 
framework to deal with the infrared sector of thermal 
theories. 



B. Counterterms 

A possible way to implement calculations in this ef- 
fective theory is to use cut-offs Jl5| ] : loops with effective 
propagators and vertices are evaluated with an upper cut- 
off A (gT < A«T). If higher order corrections are not 
necessary, the cut-off dependence is usually of a smaller 
order of magnitude. If on the contrary an important de- 
pendence on A is observed, then one must include higher 
order diagrams correcting this loop, in which A must be 
used as a lower cut-off to prevent double counting. 
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FIG. 1. Cancellation of cut-off dependence illustrated on the 
example of Baier et al. This figure also shows the 1-loop and 2-loop 
contributions when one does not use a cut-off (diagrams in boldface 
boxes). 
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One should then check that the sum of both terms is 
A independent at the desired order. This procedure has 
been used also in |l6) to calculate the hard real photon 
production rate. The results are illustrated graphically 
on figure [j], which represents numerical evaluations of 
one- and two-loop contributions to the photon polariza- 
tion tensor. 

This procedure can also be interpreted more formally 
in terms of an effective lagrangian and counterterms. The 
proper way to perform the hard thermal loop resumma- 
tion is to write an effective lagrangian £ c g which is the 
sum of the bare lagrangian C and of the term £ HTL con- 
taining all the hard thermal corrections jl7],[l8]]. There- 
fore, it appears clearly that one needs to subtract coun- 
terterms in order to keep the overall lagrangian unmodi- 
fied. Basically, this procedure amounts to write 



-£ct = ^HTL ■ 



(1) 

(2) 
(3) 



That way, the HTL resummation is just a mere reorder- 
ing of the perturbative expansion, and the exact theory 
remains the same. Counterterms have been used in |l9| ] 
for the two-loop plasmon mass. 

The subtraction of counterterms and the use of cut- 
offs are equivalent. This is illustrated on figure [I], where 
we also represented the 1-loop and 2-loop contributions 
when the cut-off is removed, in the boldface boxes. Ob- 
viously, their sum is much higher than the value of the 
plateau in the cut-off approach. This means that this 
sum includes some contributions twice, and that a coun- 
terterm is mandatory. 

Counterterms seem more systematic since they come 
from a lagrangian formulation, while it can be tricky 
to use properly cut-offs in a situation with overlapping 
loops. In the following, I follow the point of view based 
on counterterms. 



C. Photon production rate 

Let me now give the precise relationship between the 
photon production rate and the thermal Green's func- 
tions. The number of real photons emitted per unit time 
and per unit volume of the plasma is related to the re- 
tarded photon self energy via pO|,pi| 



dN 
dtcPx 



d 3 q 



2n B (g )ImII 



S(qo,q), (4) 



where n B (q D ) = l/(exp(q /T) — 1) is the Bose-Einstein 
statistical factor. For photons of invariant mass Q 2 > 
subsequently decaying in a lepton pair, we have instead 



dN 
dtd 3 x 



dq d 3 q a 
12tt 4 Q2 



Basically, the second formula differs from the first one 
by the allowed phase space, by the propagator l/Q 2 of 
a massive photon, and by the coupling constant a which 
appears in the decay into a lepton pairQ The imaginary 
part of the retarded self energy appearing in the right 
hand side can be expressed as a sum of cuts through the 
corresponding diagrams 22 24 . 



III. TWO-LOOP CONTRIBUTIONS 

A. Involved topologies 

Bremsstrahlung appears in the two-loop diagrams rep- 
resented on figure ^|. 




FIG. 2. Two-loop contributions involving bremsstrahlung pro- 
cesses. A black dot denotes an effective propagator or vertex. 
Crosses are HTL counterterms. 

According to the discussion made earlier, they should 
be accompanied by the 1-loop diagrams involving the 
proper HTL counterterms. In fact, these two- loop di- 
agrams contain much more than bremsstrahlung. For in- 
stance, if we consider a cut that goes through the gluon 
propagator, we can have bremsstrahlung if the gluon is 
space-like (L 2 < 0) and corrections to Compton effect if 
L 2 > 0. Concerning the counterterms, since their role is 
to avoid double counting, it is quite natural to associate 
them with the contributions coming from the L 2 > re- 
gion of phase space. Indeed, the inner structure of the 
counterterm involves a bare gluon which contributes only 
in the time-like region once cut. Another way to justify 
this is to say that the counterterms should be associ- 
ated with terms that provide corrections to a process 
appearing already at a lower order. There is no point in 
associating counterterms with a process that appears in 
this order for the first time. Therefore, counterterms will 



n B (g )Imn ^(q ,q) 



(5) 



1 These formulas are valid only at first order in the QED 
constant a since they neglect any reinteraction of the emit- 
ted photon on its way out of the plasma. On the contrary, 
they include corrections to all orders in the strong coupling 
constant. 



2 



be calculated only if it appears that the L 2 > region 
contributes at leading order. 

Phase space considerations tend to make us favor a 
hard quark loop in order to enlarge its size as much as 
possible. Therefore, we can drop at leading order hard 
thermal corrections involving the propagator of such a 
hard quark, and obtain a simplified version of the previ- 
ous diagrams, which is represented on figure 





FIG. 3. Simplified two-loop contributions. The circled vertices 
correspond to the framework of the cutting rules. 

More precisely, we will use the hard limit of the ef- 
fective quark propagator, retaining that way an asymp- 
totic thermal mass. This mass will be necessary later 
in order to regularize collinear singularities.^ Other au- 
thors [^,^6) have found that the width obtained by an 
extra resummation including a self energy with an imag- 
inary part in the time like region may compete with this 
asymptotic mass as a collinear regulator under certain 
circumstances. But since our goal was to follow strictly 
the HTL approach, we did not consider this extension 
here. 



B. Matrix element 

We have checked that the two cuts represented on fig- 
ure |§| form a gauge invariant set of terms. One should 
also add a second symmetrical cut to the vertex, and a 
second self-energy diagram in which the loop correction 
is on the other quark line. We can as well take them into 
account by just multiplying the result by an extra factor 
2, since they provide contributions equal to those of the 
cuts already represented. Applying the cutting rules for 
the retarded-advanced formalism p4[ gives immediately 
for the vertex contribution 



RA 

Imn /0 o ,g)|„ 



NC 



d 4 P f d A L 



,(L)-A(L) 



2 ~ a J (2tt) 4 7 (2tt) 4 
S R (P) - S A (P)l \S R (R + L) — S A (R + L) 
x {n F (r a ) - n F (p )) (n B (l a ) + n F (r a + l a )) 
x S{R)S{P + L)P T ' L {L) Traced, > 



(6) 



2 Contrary to we don't have to introduce such a mass by 
hand in the theory, since it comes naturally from the diagrams 
we are considering. 



where, we denote the fermion propagator: 



S • {P) = fS ' (P) with P = (p , vV+M^p) (7) 

S ^ " tPT" = P 2 - M 2 ± iv e ' ^ 
P ± ip a e r lvI oc ^ 'Pot 

and the effective gluon propagator in linear covariant 
gauges: 

-D pa (L) 

^ P T (L)A R T A (L) + P L (L)A R L A (L) +£VWL 4 (9) 



L 2 -U T JL) 



U T (L) = 3m 2 



x 
~2 



n . a 
x(l - X 2 ) 



In 



x 



U L (L) = '3m 2 (l - x 2 ) 



In 



x- 1 
x + 1 



X 



1 



(10) 
(11) 



with P. 
tors 



the usual transverse and longitudinal projec- 
27|H], ee g 2 C F T 2 /A [| the asymptotic ther- 
mal mass of the quark, m 2 = g 2 T 2 [N + N F /2]/9 the soft 
gluon thermal mass, and x = l /l. In this formula, e 
denote the electric charge of the quark and therefore de- 
pends on its flavor. Likewise, we obtain for the second 
diagram: 



RA 

imn /(<? ,<z)| a 
NC F 2 2 

-e g 



(If 

d A p 



d L r 



XL) 



2 ' a J (2tt) 4 J (2tt) 4 
S R (P) - S A (P)l \S R (R + L)- S A (R + L] 
x ( u f ( r o) ~ n F (p )) (n B (l ) + n F (r a + l Q )) 



(S(R)) 2 P pa (L) Trace pc 



I self 



(12) 



For propagators S(R) without any R or A superscript, 
this prescription is irrelevant because of the position of 
the cuts. It means that only the principal part of these 
propagators contributes. 

Taking into account the identity L p P pa (L) = 0, as 
well as the fact that certain terms will be killed later by 
the S(.) distributions, we can give the following expres- 
sions for the Dirac's traces 

Trace pCT | Bclf w -4 [4R 2 Q P R' J - AQ 2 R p R a 

-g pa (R 2 (R 2 -Q 2 ) + 2R 2 Q ■ L-2Q 2 R- L)] , (13) 

Trace p<T | vcrtox w -4 [2R 2 P p Q a - 2(FTL) 2 R p Q' t 
+2L 2 (R p R a + P p P a ) - AQ 2 R p P a 
+g pa \-L 2 (R 2 + (FTL) 2 - Q 2 - L 2 ))} . 



(14) 
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C. Phase space 

Let me now study the restrictions that the i5(.) func- 
tions impose on the available phase space. Since we may 
encounter collinear divergences, we need to keep carefully 
the asymptotic mass M m in this study. From the first 
cut propagator S a (P)-S A (P) = 2i:e{ Po )5{P 2 -M 2 ^), we 
obtain p Q — ±uj p , where I denote iv p = ^/(p 2 + M^), as 
well as r Q — p + q - The second delta function constrains 
the angle 9' between the 3- vectors r and I, by 

c OBe >= R2 - M ~+? r ° l ° + L2 . (15) 
2rl 

At that point, we must impose that this cosine be be- 
tween — 1 and +1, which amounts to the following set of 
inequalities 

(r - r + l + l)(r +r + l a -l)> 

{r -r + l - l)(r + r + l + I) < . (16) 

A convenient way is to see these inequalities as con- 
straints in the (l a , I) plane in which they reduce the in- 
tegration domain. The effect of this reduction is rep- 
resented on figure ^. On this figure, regions shaded in 
dark gray are excluded by the inequalities Regions 
shaded in light gray are allowed and involve a time-like 
cut gluon propagator, which means that the correspond- 
ing processes correspond to Compton or annihilation pro- 
cesses. Regions in white are also allowed and correspond 
to the exchange of a space-like gluon. These regions will 
turn out to be dominant. 



Po = -% y> 


I 


© 







\ I 










■m g -r„-(o r 






+ m g 


-r„ + a r l„ 







© 




I 




Po = +% 






III 






■ r - ra. 






■r + a r 


k 



FIG. 4. Allowed domains in the (l ,i) plane for p = ±u> p . 
The area shaded in dark gray is excluded by the delta function 
constraints. The areas shaded in light gray are above the light-cone 
(dotted lines). The light curves are the transverse and longitudinal 
mass shells of the thermalized gluon. The vertical dotted line is the 
separation between e(p )e(r + l a ) = +1 and e(p )e(r + l ) = —1. 

The regions on which we focus mostly in the following 
have been labeled by I, II and III. Looking at the signs 
of the quark momenta in each region, we can see that 
region I and III involve bremsstrahlung processes, while 
a new process appears in region II. This unexpected pro- 
cess can be described as a qq annihilation in which one 
of the fermions undergoes a scattering through a gluon 
exchange (see the figure ||) . 




FIG. 5. Physical processes included in the regions where L 2 < 0. 

By the change of variables P — > —R — L, we can map 
the region I on the region III. As a consequence, the 
bremsstrahlung photon production can be obtained by 
looking only at the region III, and multiplying its contri- 
bution by an overall factor 2. 

D. Collinear singularities 

When both the external photon and the internal quark 
are massless, the diagrams of figure ^ display a collinear 
singularity when the 3-momentum of the quark is paral- 
lel to that of the emitted photon. Moreover, contrary to 
the collinear divergences encountered in the 1-loop con- 
tributions, which are at most logarithmic, we have here 
two denominators that can vanish almost simultaneously. 
This is of course obvious for the second diagram of figure 
^ in which we have a double pole. But even for the first 
one, the poles are distinct but very close. Indeed, R 2 is 
vanishing when p is parallel to q, while (P + L) 2 van- 
ishes when r + 1 is parallel to q. As a consequence, when 
/ is soft, these two conditions of collinearity are satis- 
fied almost simultaneously, independently of the order of 
magnitude of the momentum q. Therefore, the two de- 
nominators have very close poles, and their association 
mostly behaves like a double pole. 

The collinear singularity of these diagrams is there- 
fore quasi linear, instead of logarithmic. When the quark 
mass is taken into account, these divergences are regu- 
larized, but a consequence of their linear nature is the 
generation of powers oil/g. The strength of these diver- 
gences makes dominant the terms that have this configu- 
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ration of quasi double poles. It is rather straightforward 
to show that this cannot occur when the gluon is time- 
like. Indeed, because of the delta function S(P 2 — M^), 
the denominator (P + L) 2 — can approach zero only 
if L 2 < for soft L.^ Therefore, we consider only the 
regions where L 2 < 0. It is easy to check that in the 
collinear limit, for soft L, the boundaries of region III 
are just —l<l a < I, i-e. the inequality | cos#'| < 1 does 
not add any extra constraint beyond L 2 < 0. 



E. Soft photon production 

Let me start with soft photons. I consider here photons 
having a small invariant mass Q 2 -C q 2 -C T 2 ^ In this 
situation, it is easy to verify that the contribution of re- 
gion II is subdominant, since the area of the phase space 
available to this process in the (Z D , I) plane is proportional 
to q . Moreover, for this region, the gluon momentum L 
is necessarily hard, which implies that the two conditions 
of collinearity discussed above cannot be satisfied simul- 
taneously. Therefore, bremsstrahlung is the only domi- 
nant process for the production of low invariant mass soft 
photons. 

Using only the variables u = 1 — cos 9 where is the 
angle between p and q, p, I and x — l /l, we can write 



R l - Mi 



2pq 



u + 



2p 2 



(17) 



where M 2 S = + Q 2 p 2 /q 2 , and 



2^{2 V q\ 



(P + L) 



ML 



i'ti 

~2p^ 



L 2 y _ L 2 M| ft 

Ip 1 J ~p T p' 2 



1/2 ' 



(18) 



where we performed the integration over the azimuthal 
angle 4> between q and I at this stage since this denomi- 
nator is the only place where (f> appears at the dominant 
order. Among all the terms present in the amplitude, 
the only ones to be sensitive to the collinear singular- 
ity described in the previous section are those where two 
of these denominators are present. Other terms have a 
numerator that cancels one of the denominators, so that 



3 The configuration with L 2 = and hard L has in fact been 
studied by jl||o| (there is no contribution with L 2 > in 
Jl6|]3o| ] since they used bare gluon propagators). 

4 The previous discussion on collinear singularities still holds 
for photons which are slightly virtual. The inequality Q 2 
q 2 is in fact the condition on the invariant mass in order to 
have this collinear enhancement. 



the collinear singularity is only logarithmic. On the basis 
of these considerations, it is easy to check that only the 
following term in L 2 



8L 2 ^ 



R p R a + P p P° 



(19) 



R 2 {P + L) 2 
is dominant. 

It is then straightforward to perform the angular in- 
tegration over u = 1, to obtain the following expression 
of the imaginary part of the polarization tensor of the 
photon: 

p2 9 2 N c C F T\ 7 



where we denoted 



3tt 2 



qo 



(20) 



j =1 

T o 



+ OG 



dvv 



(e v + If 



+oo 



with 



P_ 

T 



dw 



w 



^Jw/(w + 4)tanh~ 1 y/w/(w + 4) 
( W + R^ L (x)) 2 + (I T ^ L (x)) 2 



(21) 



M 2 

err 



I L (x) = 



ImIL 



(x) 



M 2 



R T Ax) 



ReII T L (x) 
M 2 , 



(22) 



A first noticeable feature of the result of Eq. (20) is 
its order of magnitude. Indeed, it is 1/g larger than 1- 
loop contributions calculated in This is due to the 

strong collinear singularities described earlier. 

The quantities J T L are coefficients quantifying the re- 
spective contributions of the transverse and longitudinal 
exchanged gluons. If we plot these quantities as a func- 
tion of the parameter Q 2 /q 2 , we obtain the curves of fig- 
ure @. One can see that these coefficients are decreasing 
very fast if Q 2 increases. This is of course due to the fact 
that the invariant mass Q 2 attenuates the effect of the 
collinear singularities. What happens around Q 2 jq 2 ~ 1 
will be presented in detail by H. Zaraket elsewhere in 
these proceedings M]. 
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FIG. 6. Transverse and longitudinal contributions as a function 
of Q 2 /q 2 , for 3 colors and 3 light flavors. The value taken for the 
coupling constant is g = 0.44. 

Another feature of the result ( pp| ) is that it is totally 
free of any infrared divergence when the gluon momen- 
tum is going to zero, while one could have expected prob- 
lems to occur with the transverse gluons (see the prob- 
lems encountered in the calculation of the damping rate 
of a moving quark, for instance). If we look formally at 
the limit — * of the quantities J T L when Q 2 = 0, 
we find 



lim J L ~ ln(mg/Moo) 

Mooting 

lim J T ~ liifwig/Moo) 2 



(23) 
(24) 



This result is interpreted as follows: the common power 
of the logarithm is due to the collinear singularity (com- 
mon to both the transverse and longitudinal gluon con- 
tributions) which reappears when its regulator is 
suppressed. The reason why only a logarithm appears in 
this limit despite the fact that the singularity has been 
described earlier as quasi linear is precisely related to the 
fact that we don't have exactly a double pole but two very 
close poles: the small separation between the two poles, 
which is responsible for extra powers of l/g, does not go 
to zero when — > 0. In the transverse gluon contri- 
bution, there is an extra power of this logarithm which 
is interpreted as a remnant of an infrared singularity. To 
support this interpretation, we can introduce by hand a 
magnetic mass m mag in the effective gluon propagator, 
and look at 



M- 



lim 

,«mma I «"'( 



J T ~ ln(m g /M 00 )ln(m g /m ma g) . (25) 



This limit indicates the infrared nature of the extra log- 
arithm in the transverse gluon contribution. Equation 
( p4| ) then tells us that the relevant regulator of this poten- 
tial infrared singularity is the quark thermal mass when 
Moo 3> m mag , which can be understood in terms of phase 
space constraints. Indeed, if L — > while > 0, 

the delta functions S(R 2 - M^) and 5((P + L) 2 - M^) 
become incompatible with bremsstrahlung, which means 
that Mqo prevents the infrared singularity from having a 
support. 



F. Hard photon production 

The calculation of 2-loop contributions to hard real 
photon production follows closely what has been done 
before for soft photons, at least for Q 2 = 0, a condition 
to which I will adhere in the following. Indeed, most of 
the aspects of the calculations were based on the fact that 
Q 2 <C q 2 . Simply, since q cannot be neglected in front of 
p, one must treat more carefully the angles between the 
various 3-vectors of the problem. If we still denote by 8 



the angle between p and q, the angle 8" between q and 
I is now given by 

cos 8" = - (cos 8 cos 8' + sin 8 sin 8' cos (f>) H — cos 8' , 
r r 

(26) 

where 8' and 4> are defined in the same way as before. 
Then, it is easy to see that (17) is unchanged (except for 
the fact that M 2 S becomes simply Ml when Q 2 = 0), 
while ([!§) becomes 



2tt 



2n{p + q){2qp 2 



(P + L) 2 - Ml 



ML 



L 2 
2pZ 



L 2 ML 

— 2 T~ 

p 2 p 2 



1/2 



(27) 



Again, the only dominant term in the matrix element is 
given by (|l^), which makes easy the evaluation of the 
bremsstrahlung contribution to hard real photons. We 
now obtain 



Imn M /(g , g ) 



e 2 g 2 NC F T f T 

A „ ( J T + J L ) 



7T 



X / dp [p 2 + (p + q ) 2 } [n F (p) - n F (p + q )] , (28) 



with now 



dx 



J T , — I |4,L ( x 



x / dw 



^Jw/w + 4 tanh 1 ^w/w + 4 
(w + R TtL (x)) 2 + (T T:L (x)) 2 



(29) 



As one can see, the structure of the bremsstrahlung con- 
tribution to the production of hard real photons is quite 
similar to the case of soft photons. In the above expres- 
sions, the coefficients J T L have been defined to match 
exactly those which have been defined previously, when 
Q 2 = 0. The only difference is in the integration over 
p which is now more complicated because we cannot ne- 
glect q . The q Q dependence that comes from this inte- 
gral is non trivial now, but some limits are simple. In the 
limit q <C T, we recover of course the result given for 
the production of soft real photons. A simple result can 
also be obtained in the opposite limit q Q T, for which 
a simple expansion can be given for the integral over p, 
which leads us to 



imn A V(^g) 



eVln(2)ATC 



T 2 (J T +J L ). (30) 



As said before, the contribution of region II is neglige- 
able only when q Q is soft because of a small phase space, 



G 



so that we must now take it into account. Its evalua- 
tion in the asymptotic region q a 3> T is very easy. Since 



Po 



in this region, we have instead of (11 



R 2 - Mi 



-2pq 



M 2 
2p 2 



(31) 



where I denote v = 1 + cos 8. We can also check that 
( p7j ) remains valid, except for the changes u — ► v and 
p + q —> q — p. The enhancement due to collinear singu- 
larities in the terms that contains these two denominators 
is unchanged near v ~ g 2 . Concerning the boundaries of 
the integration domain, we can note that the statistical 
weight n F (p ) — n F (p + q ) is equal to 1 for — q a < p Q < 
and equals everywhere else (excepted in a small region 
of width T where the transition from 1 to takes place, 
but this provides corrections of relative order T/q a <C 1). 
For the I and x variables, the limits are — 1 < x < 1 and 
< I < I* where the upper bound I* is always hard so 
that it can be taken to infinity without changing the re- 
sult. Therefore, the integration over x and I gives again 
the same factor J T + J L as in (|30| ) . Since the integration 
over p is replaced now by 



dp(p 2 + (q Q - P?) 



3 



the contribution of region II for q a T is 
e 2 g 2 NC F 



Imn /(<Zo,<7) 



3tt 4 



-q T(J T + J L 



(32) 



(33) 



It appears that the process of this region, behaving as 
q a T when q a 3> T, instead of T 2 for bremsstrahlung, is 
the dominant one for the production of extremely hard 
photons. 



IV. COMPARISON WITH EXISTING RESULTS 

A. Other thermal field theory approaches 

The production rate of soft real photons has already 
been evaluated at 1-loop in thermal field theory by 
[0,1). The typical order of magnitude they found for 
the imaginary part of the photon polarization tensor is 
e 2 g 4 \n(l/g)T 3 /q , which means that the bremsstrahlung 
contribution found in ( ^p|) is larger by a factor l/g 2 . 
Bremsstrahlung seems therefore to be the dominant pro- 
cess for the production of soft quasi-real photons. This 
extra l/g 2 factor can be traced back into the strong 
collinear singularity exhibited by bremsstrahlung pro- 
cesses. This effect can be seen as a breakdown of the 
effective expansion based on the resummation of hard 
thermal loops. Indeed, the powers of l/g generated by 
these collinear divergences can change the order in g of 
a contribution, so that the usual connection between the 
number of loops and the order in g is lost. 



Concerning hard real photons, |16|,|30| calculated in 
thermal field theory the contribution of processes involv- 
ing a time-like gluon, like Compton processes. Their cal- 
culation involves both 1-loop and 2-loop diagrams, and 
leads to 



imn 



AR 



cq Q 



(34) 



for q 3> T, where c is a numerical constant equal to 0.23. 
We are now in a position to compare the contributions 
given by bremsstrahlung, the qq annihilation appearing 
in region II, and the Compton processes calculated by 
Baier et al. For this comparison, we take N = 3 colors 
and N F = 2 light flavors, and a coupling constant g = 
0.5.^| At Q 2 — 0, the quantities J T L depend only on TV 
and N F , and are J T w 4.45 and J L w 4.26 with our choice 
of parameters. The evolution with the photon energy 
of these three contributions has been plotted on figure 
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FIG. 7. Comparison of various contributions to ImIl f ' (1 (Q) for 
a hard real photon. 

We see that bremsstrahlung is dominating for the 
smaller values of q , while the region II becomes dom- 
inant for hard enough q Q . For intermediate energies, the 
three contributions are comparable. The sum of the three 
contributions is significantly above the Compton contri- 
bution considered alone. 



B. Connection with plasmon frequency calculations 

In his talk p2| , Flechsig presented results about the 
QCD plasmon frequency near the light-cone. This calcu- 
lation is also plagued by very strong collinear singulari- 
ties, very similar to those encountered in the real photon 
rate. A first sight, there is however a difference because 



5 In the pcrturbative regime where g — > and for photon 
energies around T, Compton processes are actually dominant 
because of an extra factor of order ln(l/g). 
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the collinear singularities in the plasmon frequency arise 
already at one loop, and come from ggg and gggg effec- 
tive vertices. 

In order to make the connection between the present 
work and Flechsig's approach, it is useful to note that at 
the level of the topologies, we have the following identity 




(35) 

between our diagrams and some 1-loop tadpole diagram. 
Indeed, the tadpole contains both the self-energy and the 
vertex topologies since there are several terms contribut- 
ing to the "fjgg effective vertex, obtained by exchanging 
a photon and a gluon. Of course, the tadpole vanishes 
when one takes the trace over photon Lorentz indices, 
so that this equation does not reflect an exact algebraic 
identity between both sides (one would need to calculate 
the tadpole beyond the HTL approximation in order to 
be able to identify the two sides) . Despite this limitation, 
one can note that both sides have the same denominators 
and therefore the same singularities. As a consequence, 
we can interpret our singularities as being collinear singu- 
larities encountered in hard thermal vertices when some 
external legs are on the light-cone || and conclude that 
the singularities encountered in our calculation are ba- 
sically the same as those encountered by Flechsig. The 
only difference is that the HTL approximation makes the 
tadpole vanish when one traces it over the photon Lorentz 
indices, so that we have to go to two loops. 

This remark suggests also a defect of the HTL ap- 
proximation. Indeed, in the calculation of the tadpole 
of Eq. d35| ) the HTL approximation gives a zero coef- 
ficient to otherwise extremely singular (and large after 
regularization) terms. Therefore, it seems that the HTL 
approximation is not a good one in the presence of these 
strong collinear singularities, because one would need to 
go beyond this approximation to retain the relevant dom- 
inant terms in the tadpole. 



bremsstrahlung contribution as given by 
field theory in the following form 



dN 



dtdx 



dq 



brcmss 

R+L 



d 4 P 



{2ir) 3 2q J (2tt) 4 

> ^> h — 



in thermal 



d 4 L 



d 4 K 



(2tt) 4 J (2tt) 



I 



xn F (r a + l )n F (k )[l - n F (p )][l - n F (k + l )] 

x2nS(P 2 - Ml) 2nS((R + L) 2 - ) 

x2nS(K 2 - Ml) 2ttS((K + L) 2 - Af* ) , (36) 

where R = P + Q. This kind of formula looks like the 
starting point of the semi-classical approach, but it differs 
from it by some essential aspects. The most important 
one is that it includes the contributions of the transverse 
exchanged gluons, which we have shown to be quite im- 
portant since J T is as large as J L . Neglecting it does 
not appear to be a good approximation. Moreover, this 
formula differs from the semi-classical ones by the de- 
tails of how the Debye screening is incorporated in the 
gluon propagator. Indeed, the semi-classical approach 
just includes some constant Debye mass as a regulator, 
while in thermal field theory one takes into account all 
the momentum dependence of the gluon self-energy via 
the HTL resummation. Although this does not change 
significantly the physics, it may affect the numerical value 
of the production rate. 

On the other hand, these simplifications allow the 
semi-classical method to handle more easily the multiple 
scatterings, and to show that some new collective phe- 
nomenon may appear. More precisely, multiple coherent 
scatterings may conspire in order to give a smaller photon 
production rate, a phenomenon known as the Landau- 
Pomeranchuck-Migdal effect (3^]34|] . In order to extract 
this effect from thermal field theory, one should look at 
diagrams with two and more exchanged gluons. 



V. CONCLUSIONS AND PERSPECTIVES 



C. Semi-classical approach 

Photon production by a plasma has also been stud- 
ied by semi-classical methods in (iJ-Q]. Basically, in 
this approach, one studies the radiative energy loss of 
a fast quark undergoing multiple scatterings in a dense 
medium. The scattering centers are assumed to be static, 
which implies that the interaction is purely electric. This 
amounts to take only the longitudinal exchanged gluons 
into account. As a consequence, this interaction can be 
regularized by a Debye mass, and one does not have 
to worry about infrared divergences. In order to ex- 
hibit the analogy between this semi-classical approach 
and the thermal field theory, it is possible to write the 



In these proceedings, I have shown that 2-loop pro- 
cesses involving the exchange of a space-like gluon play a 
dominant role in the real photon production by a plasma. 
In particular, for the production of soft real photons, 
bremsstrahlung dominates over 1-loop contributions by 
a 1/g 2 factor. Therefore, this study reconciliates ther- 
mal field theory with the semi-classical approach. 

From a more formal point of view, this phenomenon 
is due to very strong collinear singularities in these pro- 
cesses. Although regularized by a thermal mass of order 
gT, an enhancement by a 1/g 2 factor appears in the re- 
sult as a remnant of these potential divergences. The pos- 
sibility of such an enhancement can make certain contri- 
butions much larger than what one would have expected 
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on the basis of their number of loops, therefore resulting 
in some kind of breakdown of the relationship between 
the order of magnitude and the number of loops. It would 
be very interesting to know if even stronger singularities 
can appear in higher order diagrams. 

There are also other motivations to look at higher or- 
der topologies. Indeed, the LPM effect found in the semi- 
classical approach seems to indicate that equally impor- 
tant contributions come from the diagrams correspond- 
ing to photon emission induced by multiple scatterings. 
Moreover, the study of the effect of the thermal width 
made by [^5 26 seems also to indicate contributions from 
higher order diagrams. 



[27] R.D. Pisarski, Physica A 158, 146 (1989). 
[28] H.A. Weldon, Phys. Rev. D 26, 1394 (1982). 
[29] N.P. Landsman, Ch.G. van Weert, Phys. Rep. 145, 141 
(1987). 

[30] J.I. Kapusta, P. Lichard, D. Seibert, Phys. Rev. D 44, 

2774 (1991). 
[31] H. Zaraket, These proceedings. 
[32] F. Flechsig, These proceedings. 

[33] L.D. Landau, I.Ya. Pomeranchuk, Dokl. Akad. Nauk. 

SSR 92, 535 (1953). 
[34] A.B. Migdal, Phys. Rev. 103, 1811 (1956). 



R. Baier, S. Peigne, D. Schiff, Z. Phys. C 62, 337 (1994). 
P. Aurenche, T. Bec herrawy, E. Pet it girard, Preprint 



ENSLAPP-A-452/93, hep-ph/940332C 



F. Flechsig, A.K. Rebhan, Nucl. Phys. B 464, 279 (1996). 
J. Cleymans, V.V. Goloviznin, K. Redlich, Phys. Rev. D 
47, 989 (1993). 

J. Cleymans, V.V. Goloviznin, K. Redlich, Z. Phys. C 
59, 495 (1993). 

R. Baier, Y.L. Dokshitzer, A.H. Mueller, S. Peigne, D. 

Schiff, Nucl. Phys. B 478, 577 (1996). 

R. Baier, Y.L. Dokshitzer, A.H. Mueller, S. Peigne, D. 

Schiff, Nucl. Phys. B 483, 291 (1997). 

P. Aurenche, F. Gelis, R. Kobes, E. Petitgirard, Phys. 

Rev. D 54, 5274 (1996). 

P. Aurenche, F. Gelis, R. Kobes, E. Petitgirard, Z. Phys. 
C 75, 315 (1997). 

P. Aurenche, F. Gelis, R. Kobes, H. Zaraket, Phys. Rev 
D 58, 085003 (1998). 

E. Braaten, R.D. Pisarski, Nucl. Phys. B 337, 569 
(1990). 

J. Frenkel, J.C. Taylor, Nucl. Phys. B 334, 199 (1990). 
P. Aurenche, T. Becherrawy, Nucl. Phys. B 379, 259 
(1992). 

MA. van Eijck, R. Kobes, Ch.G. van Weert, Phys. Rev. 
D 50, 4097 (1994). 

E. Braaten, M. Thoma, Phys. Rev. D 44, 1298 (1991). 
R. Baier, H. Nakkagawa, A. Niegawa, K. Redlich, Z. 
Phys. C 53, 433 (1992). 

E. Braaten, R.D. Pisarski, Phys. Rev. D 45, 1827 (1992). 
J. Frenkel, J.C. Taylor, Nucl. Phys. B 374, 156 (1992). 
H. Schultz, Nuc. Phys. B. 413, 253 (1994). 

H.A. Weldon, Phys. Rev. D 28, 2007 (1983). 

C. Gale, J.I. Kapusta, Nucl. Phys. B 357, 65 (1991). 

R.L. Kobes, G.W. Semenoff, Nucl. Phys. B 260, 714 

(1985) . 

R.L. Kobes, G.W. Semenoff, Nucl. Phys. B 272, 329 

(1986) . 

F. Gelis, Nucl. Phys. B 508, 483 (1997). 
A. Niegawa, Phys. Rev. D 55, 4997 (1997). 
A. Niegawa, Phys. Rev. D 56, 1073 (1997). 



9 



